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Abstract 



We consider higher-order QCD corrections to the production of colourless high- 
mass systems (lepton pairs, vector bosons, Higgs bosons, . . . ) in hadron collisions. 
We propose a new formulation of the subtraction method to numerically compute 
arbitrary infrared-safe observables for this class of processes. To cancel the infrared 
divergences, we exploit the universal behaviour of the associated transverse-momentum 
(qx) distributions in the small-gr region. The method is illustrated in general terms 
up to the next-to-next-to-leading order (NNLO) in QCD perturbation theory. As a 
first explicit application, we study Higgs boson production through gluon fusion. Our 
calculation is implemented in a parton level Monte Carlo program that includes the 
decay of the Higgs boson in two photons. We present selected numerical results at the 
LHC. 



March 2007 



The dynamics of scattering processes at high-momentum scales, Q, is well described by pertur- 
bative QCD. Thanks to asymptotic freedom, the cross section for sufficiently inclusive reactions 
can be computed as a series expansion in the QCD coupling as(Q 2 ). Until few years ago, the 
standard for such calculations was next-to-leading order (NLO) accuracy. Next-to-next-to-leading 
order (NNLO) results were known only for few highly-inclusive reactions (see e.g. Refs. [HEiE]). 

The extension from NLO to NNLO accuracy is desirable to improve the QCD predictions and 
to better assess their uncertainties. In particular, this extension is certainly important in two cases: 
in those processes whose NLO corrections are comparable to the leading order (LO) contribution; 
in those 'benchmark' processes that are measured with high experimental precision. Such a task, 
however, implies finding methods and techniques to practically achieve the cancellation of infrared 
(IR) divergences that appear at intermediate steps of the calculations. 

Recently, a new general method pi], based on sector decomposition [5], has been proposed 
and applied to the NNLO QCD calculations of e + e~ — > 2 jets [Bj, Higgs [7] and vector j8] boson 
production in hadron collisions, and to the NNLO QED calculation of the electron energy spectrum 
in muon decay j9]. The calculations of Refs. [3 [8] are encoded in publicly available numerical 
programs that allow the user to compute the corresponding cross sections with arbitrary cuts on 
the momenta of the partons produced in the final state. 

The traditional approach to perform NLO computations is based on the introduction of aux- 
iliary cross sections that are obtained by approximating the QCD scattering amplitudes in the 
relevant IR (soft and collinear) limits. This strategy led to the proposal of the subtraction [10] and 
slicing [11] methods. Exploiting the universality properties of soft and collinear emission, these 
methods were later developed in the form of general algorithms [T2| [TBI [Li] . These algorithms 
make possible to perform NLO calculations in a (relatively) straightforward manner, as soon as the 
corresponding QCD amplitudes are available. In recent years, several research groups have been 
working on general NNLO extensions of the subtraction method [151 US El EES EES]- Although 
NNLO results have been obtained only in some specific processes (e + e~ — > 2 jets [201 HI] and, 
partly, e + e~ — > 3 jets [22]), in the case of lepton collisions some of these general projects are near 
to completion. 

In this letter we reconsider the problem of the extension of the subtraction method to NNLO. 
Rather than aiming at a general formulation, we limit ourselves to considering a specific, though 
important, class of processes: the production of colourless high-mass systems in hadron collisions. 
We present a formulation of the subtraction method for this class of processes, and we apply it to 
the NNLO calculation of Higgs boson production via the gluon fusion subprocess gg — > H . This 
explicit application crosschecks the results of Ref. [7], by using a completely independent method. 

We consider the inclusive hard-scattering reaction 

h + h 2 -> F(Q) + X, (1) 

where the collision of the two hadrons hi and h 2 produces the triggered final state F. The final 
state F consists of one or more colourless particles (leptons, photons, vector bosons, Higgs bosons, 
. . . ) with momenta g« and total invariant mass Q (Q 2 = <7«) 2 )- Note that, since F is colourless, 
the LO partonic subprocess is either qq annihilation, as in the case of the Drell-Yan process, or 
gg fusion, as in the case of Higgs boson production. 

At NLO, two kinds of corrections contribute: i) real corrections, where one parton recoils 
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against F; ii) one-loop virtual corrections to the LO subprocess. Both contributions are separately 
IR divergent, but the divergences cancel in the sum. At NNLO, three kinds of corrections must 
be considered: i) double real contributions, where two partons recoil against F; ii) real-virtual 
corrections, where one parton recoils against F at one-loop order; iii) two-loop virtual corrections 
to the LO subprocess. The three contributions are still separately divergent, and the calculation 
has to be organized so as to explicitly achieve the cancellation of the IR divergences. 

Our method is based on a (process- and observable- independent) generalization of the proce- 
dure used in the specific NNLO calculation of Ref. (23]. We first note that, at LO, the transverse 
momentum q T = £\ q^j of the triggered final state F is exactly zero. As a consequence, as 
long as qT 7^ 0, the (N)NLO contributions are actually given by the (N)LO contributions to the 
triggered final state F + jet(s). Thus, we can write the cross section as 

d°~(N)NLo\qT^0 = ^ a (AT)LO • (2) 

This means that, when qT 7^ 0, the IR divergences in our NNLO calculation are those in da^Q ts : 
they can be handled and cancelled by using available NLO formulations of the subtraction method. 
The only remaining singularities of NNLO type are associated to the limit qT — > 0, and we treat 
them by an additional subtraction. Our key point is that the singular behaviour of daf^t^ when 
qT —>■ is well known: it comes out in the resummation program [24J of logarithmically-enhanced 
contributions to transverse-momentum distributions. Then, to perform the additional subtraction, 
we follow the formalism used in Ref. [251 126] to combine resummed and fixed-order calculations. 

The following sketchy presentation is illustrative; the details will appear elsewhere. We use a 
shorthand notation that mimics the notation of Ref. [25] . We define the subtraction counterterrr0 

da CT = do F LO ® Z F (q T /Q) d\ T - (3) 

The function S F (gj>/<5) embodies the singular behaviour of do~ F+sets when qT — > 0. In this limit 
it can be expressed as follows in terms of g^-independent coefficients E-^w*); 
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The extension of Eq. (J21) to include the contribution at qT = is finally: 



_ jp _ - jp - pp 

(N)NLO = 't{N)NLO ® da LQ + 



j F+jets _ 7 CT 
a(7 (N)LO a(T (N)LO 



(5) 



Comparing with the right-hand side of Eq. (J2J), we have subtracted the truncation of Eq. ([3]) at 
(N)LO and added a contribution at qT = needed to obtain the correct total cross section. The 
coefficient W/^w i0 does not depend on qT and is obtained by the (N)NLO truncation of the 
perturbative function 

H F = 1 + ?*H F W + (^) 2 H F ® + ... • (6) 

7T V 7T / 

A few comments are in order. 

• The counterterm of Eq. ([3]) regularizes the singularity of da F+>ets when qT — > 0: the term 
in the square bracket on the right-hand side of Eq. is thus IR finite (or, better, inte- 
grable over qx). Note that, at NNLO, da?^ LO acts as a counterterm for the sum of the two 



^The symbol <£> understands convolutions over momentum fractions and sum over flavour indeces of the partons. 
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contributions to da F+ i ets : the double real plus real-virtual contributions. Once da F+ i ets has 
generated a weighted 'event', da?^ LO generates a corresponding counterevent with LO kine- 
matics (i.e. with qx = 0) and with weight Y< F (qT/Q), where qx is the transverse momentum 
of F in the 'event'. 

• The explicit form of the counterterm in Eq. ([3]) has some degrees of arbitrariness. The LO 
kinematics of the counterevent can be defined by absorbing in different ways the g^-recoil 
of the 'event': the only constraint is that the 'event' kinematics smoothly approaches the 
counterevent kinematics when q? — > 0. The counterterm function T, F (qT/Q) can be defined 
in different ways: the key property is that, in the small-gr limit, it must have the form 
given in Eq. (J4]). Note that the perturbative coefficients Y, F ( n '> k ) are universal they only 
depend on the type of partons (quarks or gluon) involved in the LO partonic subprocess (qq 
annihilation or gg fusion). 

• The simplicity of the LO subprocess is such that final-state partons actually appear only in 
the term da F+ ^ ets on the right-hand side of Eq. (jSJ). Therefore, arbitrary IR-safe cuts on the 
jets at (N)NLO can effectively be accounted for through a (N)LO computation. Owing to 
this feature, our NNLO extension of the subtraction formalism is observable-independent. 

• At NLO (NNLO), the physical information of the one-loop (two-loop) virtual correction to 
the LO subprocess is contained in the coefficients TiS 1 ' (H^). Once an explicit form of 
Eq. ([3]) is chosen, the hard coefficients Tl F ( n > are uniquely identified (a different choice would 
correspond to different T-C F ^). 

According to Eq. (jSJ), the NLO calculation of do F requires the knowledge of TC F ^' and the LO 
calculation of da F+ ^ ets . The general (process-independent) form of the coefficient TC F ( 1 > is basically 
known: the precise relation between TL F ^ and the IR finite part of the one-loop correction to a 
generic LO subprocess is explicitly derived in Ref. |27j . 

At NNLO, the coefficient H F ^ is also needed, together with the NLO calculation of da F+ ^ ets . 
Although the general structure^ of the coefficients / H F ^ 1 is presently unknown, we have completed 
the calculation of Ti, H ^ for Higgs boson production in the large-M top limit. Since the NLO 
corrections to gg — > H + jet(s) are available [28] in the same limit, we are able to present a first 
application of Eq. (j3J) at NNLO. We have encoded our computation in a parton level Monte Carlo 
program, in which we can implement arbitrary IR-safe cuts on the final state. 

In the following we present numerical results for Higgs boson production at the LHC. We use 
the MRST2004 parton distributions [2H], with densities and «s evaluated at each corresponding 
order (i.e., we use (n+l)-loop «s at N n LO, with n = 0, 1, 2). The renormalization and factorization 
scales are fixed to the value Hr — — Mh, where Mh is the mass of the Higgs boson. 

In Fig. [1] we consider Mh = 125 GeV, and we show the bin-integrated rapidity distribution 
of the Higgs boson at LO (dotted), NLO (dashed) and NNLO (solid). The impact of the NNLO 
corrections on the NLO result is mildly dependent on the rapidity yn when \uh \ ~3. The total 
cross section increases by about 19% when going from NLO to NNLO. 

^More precisely, the NNLO coefficients E^ 2 ' 1 ) and Y> F( ~ 2 > 2) have a non-universal contribution that, nonetheless, 
is proportional to the NLO coefficient TL Ftyl \ 

§It could be derived by extending the 0{a\) calculation of Ref. [27] to compute subleading logarithms. Work 
along these lines is under way. 
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Figure 1: Bin-integrated rapidity distribution of the Higgs boson with Mh = 125 GeV: results at 
LO (dotted), NLO (dashed) and NNLO (solid). 



When searching for the Higgs boson in the H — > WW channel, a jet veto is typically required 
to suppress the WW background from ti production. In Fig.[2]we present the rapidity distribution 
of the Higgs boson with Mh = 165 GeV. In this case we apply a veto on the jets that recoil against 
the Higgs boson. Jets are reconstructed by using the kx algorithm [301 with jet size D = 0.4 S 
each jet is required to have transverse momentum smaller than 40 GeVlll As is known (23J [7] , the 
impact of higher-order corrections is reduced when a jet veto is applied. In the present case, the 
impact of the NNLO corrections on the NLO total cross section is reduced from 20 to 5 %. 

We finally consider the Higgs boson decay in the H — > 77 channel and follow Ref. [32J to apply 
cuts on the photons. For each event, we classify the photon transverse momenta according to their 
minimum and maximum value, pxmin and Prmax- The photons are required to be in the central 
rapidity region, \q\ < 2.5, with prmm > 35 GeV and prmax > 40 GeV. We also require the photons 
to be isolated: the hadronic (partonic) transverse energy in a cone of radius R = 0.3 along the 
photon direction has to be smaller than 6 GeV. When Mh = 125 GeV, by applying these cuts the 
impact of the NNLO corrections on the NLO total cross section is reduced from 19% to 11%. 

In Fig. [3] we plot the distributions in pr-min and pTmax for the gg — > H — > 77 signal. We 
note that the shape of these distributions sizeably differs when going from LO to NLO and to 
NNLO. The origin of these perturbative instabilities is well known [33J. Since the LO spectra 
are kinematically bounded by px < M#/2, each higher-order perturbative contribution produces 
(integrable) logarithmic singularities in the vicinity of that boundary. More detailed studies are 

'in our calculation up to NLO, the kr algorithm and the cone algorithm [31] are equivalent. At NNLO, the 
kr algorithm is equivalent to the cone algorithm (with cone size R = D) without midpoint seeds, while the cone 
algorithm with midpoint seeds would lead to (slightly) different results. The cone algorithm without midpoint seeds 
would be infrared unsafe starting from N 3 LO. 

"At NNLO, a jet may consist of two partons. In this case, the transverse momentum of the jet is the vector 
sum of the transverse momenta of the two partons. 
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Figure 2: Bin-integrated rapidity distribution of the Higgs boson with Mh = 165 GeV. Final-state 
jets are required to have transverse momentum smaller than 40 GeV. 

necessary to assess the theoretical uncertainties of these fixed-order results and the relevance of 
all-order resummed calculations. A similar comment applies to the distribution of the variable 
(prmin +PTmax)/2, which is computed, for instance, in Refs. [7] IM]. 

We have illustrated an extension of the subtraction formalism to compute NNLO QCD correc- 
tions to the production of high-mass systems in hadron collisions. We have considered an explicit 
application of our method to the NNLO computation of gg —>■ H —>■ 77 at the LHC, and we have 
presented few selected results, including kinematical cuts on the final state. The computation 
parallels the one of Ref. [7], but it is performed with a completely independent method. In the 
quantitative studies that we have carried out, the two computations give results in numerical 
agreement. In our approach the calculation is directly implemented in a parton level event gener- 
ator. This feature makes it particularly suitable for practical applications to the computation of 
distributions in the form of bin histograms. We plan to release a public version of our program in 
the near future. We also plan to apply the method to other hard-scattering processes. 

Acknowledgements. We would like to thank Daniel de Florian for helpful discussions and 
comments. 
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